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Since large pr particles in high-energy hadronic or nuclear collisions come from jet fragmentation, 
jet quenching due to parton energy loss in dense matter will cause the suppression of large pr hadron 
spectra in high-energy heavy-ion collisions. Assuming an effective energy loss dE/dx for the high 
Et partons, effective jet fragmentation functions are constructed in which leading hadrons will 
be suppressed. Using such effective fragmentation functions, high pr hadron spectra and particle 
suppression factors relative to pp collisions are estimated in central high-energy nuclear collisions 
with a given range of the assumed dE/dx. It is found that the suppression factors are very sensitive 
to the value of the effective energy loss. Systematic nuclear and flavor dependence of the hadron 
spectra are also studied. 
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I. INTRODUCTION 

An ideal quark-gluon plasma (QGP) has often been 
defined as a system of weakly interacting quarks and 
gluons in both thermal and chemical equilibrium. How- 
ever, recent theoretical investigations based on pertur- 
bative QCD-inspired model show that it is increas- 
ingly difficult for the initially produced partons to evolve 
into thermal equilibrium, let alone chemical equilibrium. 
Therefore, one might have a generalized QGP simply as 
an interacting and deconfined parton system with a large 
size and long life-time. 

One can find many examples of an interacting parton 
system in collisions involving strong interaction. But so 
far none of them can be considered a QGP in terms of 
either the ideal or generalized definition. At a distance 
much smaller than the confinement scale Aqcd and nor- 
mally in the earliest time of the collision, the interac- 
tion can be described by perturbative QCD (pQCD). 
Later on, the produced partons will then combine with 
each other via non-perturbative interactions and finally 
hadronize into hadrons. Therefore, one can consider 
that there exists an interacting parton system during the 
prehadronization stage in, e.g., e"'"e~ annihilation and 
deeply inelastic e~p processes, which is, however, limited 
only to a space-time region characterized by the confine- 
ment scale Aqcd- The characteristic particle spectrum 
(in pt and rapidity) and the ratios of produced particles 
are then determined by the physics of pQCD and non- 
perturbative hadronization. In ultrarelativistic heavy- 
ion collisions, one seeks to produce a similar interacting 
parton system but at a much larger scale of the order 
of a nucleus size and for a long period of time (e.g., a 
QGP). Therefore, one should study those experimental 



observables which are unique to the large size and long 
life-time of an interacting partonic system as signals of 
a quark-gluon plasma. 

Among many proposed signals of a quark-gluon 
plasma hard probes associated with hard processes 
are especially useful because they are produced in the 
earliest stage of the collision and their abilities to probe 
the dense matter are less complicated by the hadroniza- 
tion physics. Merits of hard probes are even more appar- 
ent at high energies because those processes also domi- 
nate the underlying collision dynamics which will deter- 
mine the initial conditions of the produced partonic sys- 
tem . Study of them will then enable us to probe the 
early parton dynamics and the evolution of the quark- 
gluon plasma. 

In general, one can divide the hard probes into two 
categories: thermal emission and particle suppression by 
the medium. Particle production, like photon/dilepton 
and charm particles, from thermal emission can be con- 
sidered as thermometers of the dense medium. Their 
background comes from the direct production in the ini- 
tial collision processes. On the other hand, suppression 
of particles produced in the initial hard processes, like 
high-pT particles from jets and J/^*, can reveal evidences 
of the parton energy loss in dense matter and the decon- 
finement of the partonic system. Thermal production of 
these particles is expected to be negligible. Therefore, 
in both cases, one needs to know the initial production 
rate accurately enough. Another advantage of these hard 
probes is that the initial production rate can be calcu- 
lated via pQCD, especially if we understand the modest 
nuclear modification one would expect to happen. 

In this paper, we will discuss high-p^ particles as 
probes of the dense matter since one expect high-i?T 
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partons which produce these high-p^ particles will in- 
teract with the dense medium and lose energy. Medium- 
induced energy loss of a high-energy parton traversing 
a dense QCD medium is interesting because it depends 
sensitively on the density of the medium and thus can 
be used as a probe of the dense matter formed in ul- 
trarelativistic heavy-ion collisions. As recent studies 
demonstrated 1^-0] , it is very important to take into ac- 
count the coherent effect in the calculation of radiation 
spectrum induced by multiple scattering of a fast par- 
ton. The so-called Landau-Pomeranchuk-Migdal effect 
can lead to very interesting, and sometimes nonintuitive 
results for the energy loss of a fast parton in a QCD 
medium. Another feature of the induced energy loss is 
that it depends on the parton density of the medium via 
the final transverse momentum broadening that the par- 
ton receives during its propagation through the medium. 
One can therefore determine the parton density of the 
produced dense matter by measuring the energy loss of 
a fast parton when it propagates through the medium. 

Unlike in the QED case, where one can measure di- 
rectly the radiative energy loss of a fast electron, one 
cannot measure directly the energy loss of a fast leading 
parton in QCD. Since a parton is normally studied via a 
jet, a cluster of hadrons in the phase space, an identified 
jet can contain particles both from the fragmentation of 
the leading parton and from the radiated partons. If we 
neglect the px broadening effect, the total energy of the 
jet should not change even if the leading parton suffers 
radiative energy loss. What should be changed by the 
energy loss are the particle distributions inside the jet or 
the fragmentation functions and the jet profile. There- 
fore, one can only measure parton energy loss indirectly 
via the modification of the jet fragmentation functions 
and jet profile. For this purpose, it was recently proposed 
[^J^ that the jet quenching can be studied by measuring 
the pt distribution of charged hadrons in the opposite di- 
rection of a tagged direct photon. Since a direct photon 
in the central rapidity region {y = 0) is always accom- 
panied by a jet in the opposite transverse direction with 
roughly equal transverse energy, the pr distribution of 
charged hadrons in the opposite direction of the tagged 
direct photon is directly related to the jet fragmenta- 
tion functions with known initial energy. One can thus 
directly measure the modification of the jet fragmenta- 
tion and then determine the energy loss suffered by the 
leading parton with given initial energy. 

Similarly, single-particle spectrum can also be used to 
study the effect of parton energy loss as first proposed 
in Ref. jl^, since the suppression of large Et jets natu- 
rally leads to the suppression of large pt particles. How- 
ever, since the single-particle spectrum is a convolution 
of the jet cross section and jet fragmentation function, 
the suppression of produced particles with a given px re- 
sults from jet quenching with a range of initial transverse 
energies. Therefore, one cannot measure the modifica- 



tion of the jet fragmentation function or the energy loss 
of a jet with known initial transverse energy from the 
single-particle px spectrum as precisely as in the case 
of tagged direct photons. One clear advantage of single 
inclusive particle spectrum is the large production rate 
of moderately high px particles, while the production 
rate of large px direct photons is relatively much smaller 
at the designed luminosity of the Relativistic Heavy-Ion 
Collider (RHIC) 1^. Therefore, with much less experi- 
mental effort, one can still study qualitatively the effect 
of jet quenching and extract the average value of the 
parton energy loss from single particle spectra at high 
Px- 

In this paper, we will conduct a systematic study of 
the effects of parton energy loss on single-particle trans- 
verse momentum spectra in central A-l-A collisions in the 
framework of modified effective jet fragmentation func- 
tions. We study within this framework the dependence 
of the spectra on the effective parton energy loss. We will 
discuss the energy or px and A dependence of the energy 
loss and jet quenching. Finally, fiavor dependence of the 
spectra will be also be discussed. 



II. MODIFIED JET FRAGMENTATION 
FUNCTIONS 

Jet fragmentation functions have been studied exten- 
sively in e~^e~ , ep and pp collisions ||ll[] . These functions 
describe the particle distributions in the fractional en- 
ergy, z — Eh/Ejet, in the direction of a jet. The mea- 
sured dependence of the fragmentation functions on the 
momentum scale is shown to satisfy the QCD evolution 
equations very well. We will use the parameterizations 
of the most recent analysis in both z and for 

jet fragmentation functions D^^^{z,Q^) to describe jet 
(a) fragmentation into hadrons {h) in the vacuum. 

In principle, one should study the modification of jet 
fragmentation functions in a pcrturbative QCD calcula- 
tion in which induced radiation of a propagating parton 
in a medium and Landau-Pomeranchuk-Migdal interfer- 
ence effect can be dynamically taken into account. How- 
ever, for the purpose of our current study, we can use a 
phenomenological model to describe the modification of 
the jet fragmentation function due to an effective energy 
loss dE/dx of the parton. In this model we assume: (1) 
A quark-gluon plasma is formed with a transverse size of 
the colliding nuclei, Ra- A parton with a reduced energy 
will only hadronize outside the deconfined phase and the 
fragmentation can be described as in e+e^ collisions. (2) 
The mean-free-path of inelastic scattering for the parton 
a inside the QGP is \a which we will keep a constant 
through out this paper. The radiative energy loss per 
scattering is ta- The energy loss per unit distance is 

thus dEa/dx = Ca/Xa- 
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The probability for a parton to scatter n times within 
a distance AL is given by a Poisson distribution, 



Pa{n,AL) = 



(1) 



We also assume that the mean-free-path of a gluon is 
half that of a quark, and the energy loss dE/dx is twice 
that of a quark. (3) The emitted gluons, each carrying 
energy on the average, will also hadronize according 
to the fragmentation function with the minimum scale 
Qq = 2.0 GeV^. We will also neglect the energy fluc- 
tuation given by the radiation spectrum for the emitted 
gluons. Since the emitted gluons only produce hadrons 
with very small fractional energy, the final modified frag- 
mentation functions in the moderately large z region are 
not very sensitive to the actual radiation spectrum and 
the scale dependence of the fragmentation functions for 
the emitted gluons. 

This is definitely a simplified picture. In a more realis- 
tic scenario, one should also consider both the longitudi- 
nal and transverse expansion. Because of the expansion, 
the actual parton energy loss will change as it propagates 
through the evolving system resulting in a different total 
energy loss as recently studied in Ref. jl4j. Since we are 
mostly interested in the overall effects, we can neglect 
the details of the evolution history and concentrate on 
the modification of high px hadron spectra due to an as- 
sumed total energy loss or averaged energy loss dE/dx 
per unit distance. It might require much more elabo- 
rated study to find out the effects of the dependence of 
the energy loss on the dynamical evolution of the system. 
It is beyond the scope of this paper. 

We will consider the central rapidity region of high- 
energy heavy-ion collisions. We assume that a parton 
with initial transverse energy Et will travel in the trans- 
verse direction in a cylindrical system. With the above 
assumptions, the modified fragmentation functions for a 
parton traveling a distance AL can be approximated as. 



A./a(^, Q', AL) = ^ ^ Pain, AL)'^Dl,,{zl, Q^) 

(2) 



/ 
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+ {na)^Dl,^{z',,Qi), 



where z^ = zjiX - nea/Er), 4 = zEx/ea and C% = 
Sn=o Pain). We limit the number of inelastic scattering 
to = Er/ea by energy conservation. For large values 
of iV, the average number of scattering within a distance 
AL is approximately (ria) ~ AL/Xa- The first term 
corresponds to the fragmentation of the leading partons 
with reduced energy Ex — nea and the second term comes 
from the emitted gluons each having energy on the 
average. Detailed discussion of this modified effective 
fragmentation function and its limitations can be found 
in Ref. il. 



III. ENERGY LOSS AND SINGLE-PARTICLE Pt 
SPECTRUM 

To calculate the px distribution of particles from jet 
fragmentation in pp and central heavy-ion collision, one 
simply convolutes the fragmentation functions with the 
jet cross sections |lq]. 



dyd'^PT 



for pp and 
dyd^PT 



abcdh 



dXn 



dXbfa/piXa,Q^) 



, , ^2, ^l/cM') da , ^ 
Jb/pixb,Q ) ;7(ao 



dt 



K j d\t\ir) J2 f dxa f 

abcdh-' ''-"-i^ "'^ 
fa/AiXa,Q^,r) hi A iXb,Q^,r) 

Dh/cizc,Q'^,^L)da ^ ^ 
-{ab cd), 



cd), (3) 



dxh 



(4) 



for AA collisions, where Zr 



xrieVxa + e-y/xb)/2, 

Xbmin = XaXTe~y /i2Xa - Xxe^), Xamin = XtB^ / i2 - 

XTe~y), and xt = 2ptI \fs. The nuclear thickness func- 
tion is normalized to / (i^rtAi^) = A. The K 2 factor 
accounts for higher order corrections The parton 

distributions per nucleon in a nucleus (with atomic mass 
number A and charge number Z), 



fa/Aix, Q^, r) = Sa/Aix, r) 



j^fa/pix, Q 



il-j)fa/nix,Q') 



(5) 



are assumed to be factorizable into parton distributions 
in a nucleon fa/Ni^, Q^) and the parton shadowing fac- 
tor Sa/Aix,r) which we take the parameterization used 
in HIJING model jl7| . Neglecting the transverse expan- 
sion, the transverse distance a parton produced at (r, 0) 
will travel is AL(r, (j)) = ^ i?^ — (1 — cos^ (p) ~r cos (p. 

In principle, one should also take into account the in- 
trinsic transverse momentum and the transverse momen- 
tum broadening due to initial multiple scattering. These 
effects (so-called Cronin effects) are found very impor- 
tant to the final hadron spectra at around SPS energies 
(Ts = 20 - 50GeV) However, at RHIC energy 

which we are discussing in this paper, one can neglect 
them (about 10-30% correction) to a good approxima- 
tion. 

We will use the MRS D— / parameterization of the par- 
ton distributions |19| in a nucleon. The resultant px 



spectra of charged hadrons (tt^, K^) for pp and pp col- 
lisions are shown in Fig. ^ together with the experimen- 
tal data |^-§2| for = 63, 200, 900 and 1800 GeV. 
The calculations (dot-dashed line) from Eq. (0) with the 
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jet fragmentation functions given by Ref. [ |12yi3[ | agree 
with the experimental data remarkably well, especially 
at large pt- However, the calculations are consistently 
below the experimental data at low p-r, where we believe 
particle production from soft processes, like string frag- 
mentation of the remanent colliding hadrons, becomes 
very important. To account for particle production at 
smaller pT, we introduce a soft component to the parti- 
cle spectra in an exponential form. 



dycPpT 



(6) 



with a parameter T = 0.25 GeV/c. This exponential 
form is a reasonable fit to the data of hadron pT spectra 
of collisions at = 200 GeV below px < 2 GeV/c. 
The fit is not very good below px = 0.5 GeV and the 
parameter T should also depend on colliding energy ^/s. 
However, for a rough estimate of the spectra at low pT 
this will be enough and we will keep T a constant. 

The normalization in Eq. (^) is determined from the 
charged hadron rapidity density in the central region; 



C 



1 



hard 



27rT2 V dy 



where 



dy 



1 



dy 



J2„ d'^hard 



PT 



dyd^PT ' 



(7) 



(8) 



Table | lists the values of the charged hadron rapidity 
density and the inelastic cross sections of pp collisions 
from HIJING calculations which we will use to determine 
the normalization in Eq. at different energies. 

The total px spectrum for charged hadrons in pp col- 
lisions including both soft and hard component is then. 



dNPP 



df^Zrd 



dyd^PT dyd^pT crf^ dyd'^pr 



(9) 



which are shown in Fig. ^ as solid lines. As one can see 
it improves the agreement with data at lower transverse 
momentum. 
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FIG. 1. The charged particle pr spectra in pp and pp 
collisions. The dot-dashed lines are from jet fragmentation 
only and solid lines include also soft production parametrized 
in an ex po nential form. The experimental data are from 
Ref. 0-0. 



We now also assume that the charged multiplicity from 
soft particle production is proportional to the total num- 
ber of wounded nucleons in AA collisions which scales 
like A, while the production from hard processes is pro- 
portional to the number of binary nucleon-nucleon col- 
lisions which scales like A"^^^. At low pr both type of 
processes contribute to the particle spectrum. Therefore 
the A-scaling of the spectrum at low pr depends on the 
interplay of soft and hard processes. In HIJING model 
[ p^ with a cut-off of Eto = 2 GeV for jet production 
the low-pT spectra scale like A^-^. To take into account 
of the uncertainty due to the interplay between soft and 
hard processes, we assume hadron spectrum in central 
A A collisions is, 





(GeV) 


63 


200 


900 


1800 


dN^P/dy 


1.9 


2.4 


3.2 


4.0 


oil (mb) 


35 


44 


50 


58 





TABLE I. Charged hadron rapidity density and inelastic 
cross sections for pp collisions at different colliding energies 
from HIJING calculations 



dN 



AA 



= A°"^ 



dNPP, 

soft 



dN, 



AA 



dyd'^PT dyd^pT ' 



(10) 



1.1. 



dyd^PT 

where an = 1-0 ^ 

To calculate dN^^^Jdyd'^ Pt, we will take into account 
both the effect of nuclear shadowing on parton distribu- 
tions and the modification of the jet fragmentation func- 
tions due to parton energy loss inside a medium. From 
Eq. we see that it will be proportional to overlap func- 
tion of central AA collisions Taa{^)- In a hard sphere 
model for nuclear distribution, Ta^(O) = QA'^ /^tiR\ and 
Ra = 1.2A1/3 fm. 
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We now define an effective suppression factor, or the 
ratio, 



Raa [pi 



dNAA/dy/(PpT 
afPTAAiO)dNpp/dy/d^PT 



(11) 



between the spectrum in central AA and pp colhsions 
which is normahzed to the effective total number of bi- 
nary NN collisions in a central AA collision. If none 
of the nuclear effects (shadowing and jet quenching) are 
taken into account, this ratio should be unity at large 
transverse momentum. Shown in Fig. ^ are the results 
for central Au + Au collisions at the RHIC energy with 
dEg/dx = 1, 2 GeV/fm, and A, = 1 (sohd), 0.5 fm 
(dashed), respectively. As we have argued before, jet 
energy loss will result in the suppression of high px par- 
ticles as compared to pp collisions. Therefore, the ratio 
at large pr in Fig. || is smaller than one due to the energy 
loss suffered by the jet partons. It, however, increases 
with pt because of the constant energy loss we have as- 
sumed here. At hypothetically large pt when the total 
energy loss is negligible compared to the initial jet en- 
ergy, the ratio should approach to one. 
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dEq/dx=2GeV/fm 
Au+Au(b=0) s^'^=200 GeV 



2 4 6 8 10 12 14 16 18 20 
(GeV/c) 

FIG. 2. The suppression factor or ratio of charged par- 
ticle Pt spectrum in central Au + Au over that of pp col- 
lisions at -y/s = 200 GeV , normalized by the total bi- 
nary nucleon-nucleon collisions in central Au + Au colli- 
sions, with different values of the energy loss dEq/dx and 
the mean-free-path Xq of a quark inside the dense medium, 
the shaded area indicates the uncertainty of the effective 
j4-scaling of low-pT spectra depending on the interplay of 
soft and hard processes 

Since there is always a coronal region with an average 
length of Xq in the system where the produced parton 
jets will escape without scattering or energy loss, the 
suppression factor can never be infinitely small. For the 
same reason, the suppression factor also depends on the 



parton's mean- free-path, A^. It is thus difficult to extract 
information on both dEg/dx and Xg simultaneously from 
the measured spectra in a model independent way. 

At small Pt, particles from soft processes (or from 
hadronization of QGP) dominate. The ratio Raa{pt) 
is then very sensitive to the A-scaling behavior of the 
soft particle production. Since we assumed an effec- 
tive scaling, A"'^ with ah = 1.0 ~ 1.1, for the low- 
Pt particle production, the ratio should approach to 
A^'-fappTAAiO) = 0.149 - 0.253 at small pr for cen- 
tral Au + Au collisions at the RHIC energy, as shown 
in Fig. ^. Therefore, the shaded area in the figure (we 
only plotted for one case of energy loss) should be con- 
sidered as one of the uncertainties of the ratio at low px 
associated with the interplay of contributions from soft 
and hard processes. One presumably can determine this 
dependence from future RHIC experimental data. 

As we have stated earlier, the Cronin effect due to 
initial multiple parton scattering will introduce an un- 
certainty of 10-30% effect which can be narrowed down 
through a systematic study ofp + p and p + A collisions. 
Since the concept of parton energy loss can only be ap- 
plied to high Et jets, it will in principle only affect the 
spectra at high px where contribution from soft produc- 
tion is negligible. At smaller values of < 3 ~ 4 GeV/c 
where soft particle production become important, the 
connection between parton energy loss and the hadron 
spectra becomes unclear. In this region, the modification 
of the spectra is driven by parton and hadron thcrmaliza- 
tion. Since we approximate the spectra in this region by 
an effective A-scaling of coherent or semi-coherent par- 
ticle production, the suppression factor we show in this 
paper can only be considered as semi-quantitative. Fur- 
thermore, the spectra just above pT ^ 4: GeV/c should 
also be sensitive to the energy dependence of the energy 
loss as we will show in the next section. 

To further illustrate the effect of the parton energy 
loss in hadron spectrum we show in Fig. |3| the produc- 
tion rates of vr*^ with (solid line) and without parton en- 
ergy loss (dashed lines), together with the spectrum of 
direct photons (dot-dashed lines) at the RHIC energy. 
The upper curve for direct photons is a leading order 
calculation multiplied by K — 2 factor. The lower curve 
is the result of a next-to-leading order calculation ||2^ 
which also includes quark bremsstrahlung. Here we as- 
sumed the \ow-pT soft particle spectra scales like A^-^. 
Since we can neglect any electromagnetic interaction be- 
tween the produced photon and the QCD medium, the 
photon spectrum will not be affected by the parton en- 
ergy loss. On the other hand, jet quenching due to par- 
ton energy loss can significantly reduce tt*^ rate at large 
Pt- Therefore the change of j/tt'^ ratio at large pT can 
also be an indication of parton energy loss. One can 
consider the contribution to direct photon production 
from bremsstrahlung as quark fragmentation into a pho- 
ton, it should in principle also be affected by the quark 
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energy loss inside the dense medium. Therefore, there 
is also some uncertainties ( maximum factor of 2 if the 
K — 2 factor completely comes from bremsstrahlung 
correction) to the estimated photon spectra at lower pT 
where bremsstrahlung is more important. 
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FIG. 3. The inclusive pr distribution for tt" witli (solid) 
and without (dashed) parton energy loss as compared to that 
of direct photons (dot-dashed) in central Au + Au collisions 
at y/s — 200 GeV. dEq/dx — 1 GeV/fm and mean-free-path 
Ag = 1 fm are assumed. Contribution from soft particle pro- 
duction to the TTo spectra is assume to have a A^'^ scaling. 



IV. ENERGY AND A DEPENDENCE OF 
ENERGY LOSS 



In recent theoretical studies of parton energy loss 
it has been demonstrated that the so-called 
Landau-Pomeranchuk-Midgal (LPM) coherent effect can 
lead to interesting and sometimes nonintuitive results. 
Baier et al have systematically studied these effects in 
detail ^j^- They found that because of the modifica- 
tion of the radiation spectrum by the LPM coherence, 
the energy loss experienced by a fast parton propagat- 
ing in an infinite large medium has a nontrivial energy 
dependence, 



dE 
dx 



(fori > Lc 



(12) 



where — i, E parton's energy, /i^ the Debye screen- 
ing mass for the effective parton scattering, A parton's 
mean-free-path in the medium, and L^r = \/ \E / ji^. 
For a more energetic parton traveling through a medium 



with finite length (L < L^c), the final energy loss be- 
comes almost independent of the parton energy and can 
be related to the total transverse momentum broadening 
acquired by the parton through multiple scattering. 



dE 

dx 



— T^^Pt 



Nca, 



(13) 



where Sp^ is the transverse momentum kick per scatter- 
ing the parton acquires during the propagation. There- 
fore, the energy loss per unit distance, dE/dx, is propor- 
tional to the total length that the parton has traveled. 
Because of the unique coherence effect, the parton some- 
how knows the history of its propagation. 

These are just two extreme cases of parton energy and 
the medium length. Since it involves two unknown pa- 
rameters of the medium, it is difficult to determine which 
case is more realistic for the system of dense matter pro- 
duced in heavy-ion collisions. We will instead study the 
phenomenological consequences of these two cases in the 
final single inclusive particle spectrum at large px ■ 
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FIG. 4. The same as Fig. except that an en- 
ergy-dependent energy loss is assumed. The mean-free-path 



Shown in Fig. ^ are the calculated suppression factors 
with an energy-dependent parton energy loss, dEq/dx = 
y/W/bGeV and dEjdx = y^E/20GeY GeV/fm, respec- 
tively for central Au + Au collisions at the RHIC energy. 
Comparing to Fig. || with a constant energy loss, sup- 
pression factors are fiatter as functions of pT- This is 
understandable because the energy loss for larger Et jet 
will lose more energy in this scenario thus leading to a 
stronger suppression of high particles. As pointed out 
in Ref . |^ , the most relevant quantity in the modification 
of the fragmentation functions is the parton energy loss 
AEt relative to its original energy Et- For a constant 
energy loss AEt, the ratio AEt/Et becomes smaller 



6 



for larger Et, thus the suppression factor Raa{pt) will 
increase with px- If the energy loss AEt rises with the 
initial energy Et, then the increase will be slower. Thus 
the slope of the ratio RAAiPx) can provide us informa- 
tion about the energy dependence of the energy loss, as 
one can see from the comparison of Figs. ^ and ^. 

0.6 



_ 0.5 
6 0.4 
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dEq/dx=0.5 GeV/fm 



dEq/dx=0.2 (L/fm) GeV/ffi 



A+A(b=0) s'^^=200GeV 



0.3 
0.2 
0.1 




^1/3 

FIG. 5. The suppression factor for central A + A collisions 
at pt = 10 GeV/c, as a function of the system size, A^^^ , for 
a constant energy loss per unit distance length (dot-dashed 
line) and an energy loss which increases linearly with the 
length (solid line). The mean-free-path A, = 1 fm is used. 

To study the consequences of a parton energy loss 
dE/dx which increases with the distance L it travels, one 
can either vary the impact-parameter or the atomic mass 
of the projectile and target so as to change the size of the 
dense matter through which the leading partons have to 
propagate. Assuming a transverse size of the colliding 
nuclei which have a hard sphere nuclear distribution, one 
can estimate that the averaged distance a produced jet 
has to travel through is (L)a = 1.09A^^^, where one has 
to weight with the probability of jet production or the 
overlapping functions of A A collisions. In Fig. we plot 
the suppression factor Raa{pt) at a fixed px—lO GeV/c 
for central A + A collisions at the RHIC energy as a 
function of A^/"^, A the atomic masses of the projectile 
and target nuclei. The solid line is for an energy loss, 
dE/dx = 0.2(L/fm) GeV/fm, which is proportional to 
the total length traveled by the parton, and dot-dashed 
line is for a constant dE/dx = 0.5 GeV/fm. As the size 
of the system increases, a parton will lose more energy 
and thus will lead to increased suppression in both cases. 
For collisions of heavy nuclei ( A^^^ > 3), the energy loss 
in the first case becomes larger than the second one and 
thus leads to more suppression. However, the functional 
form of the A-dependence of the suppression factor in 
the two cases do not differ dramatically. It is therefore 



difficult to determine whether the energy loss per unit 
length is proportional to the total length simply from 
the A-dependence of the suppression factor. It must re- 
quire a model dependent phenomenological study of the 
experimental data. 



V. FLAVOR DEPENDENCE 



Because of the non-abelian coupling, gluons in QCD 
always have stronger interaction than quarks. The gluon 
density inside nucleons at small x is larger than quarks; 
gluon-gluon scattering cross section is larger than quark- 
quark; and a gluon jet produces more particles than a 
quark jet. For the same reason, a high energy gluon will 
also lose more energy than a quark propagating through 
a dense medium. Theoretical calculations all show 
that gluons lose twice as much energy as quarks. In this 
section we will discuss how to observe such difference in 
the final hadron spectrum. 

By charge and other quantum number conservation, 
fragmentation functions of a gluon jet into particle and 
anti-particle will be identical, though it produces more 
particles than a quark jet and consequently its fragmen- 
tation functions are often softer than a quark's, as has 
been measured in the three-jet events of e"'"e~ annihi- 
lation pS) . For example, equal number of protons and 
anti-protons will be produced in the gluon fragmenta- 
tion. On the other hand, an up or down quark is more 
likely to produce a leading proton than anti-proton and 
vice versa for anti-quarks. Since there will be more quark 
(up and down) jets produced than anti-quark in nuclear 
collisions, one will find more protons than anti-protons, 
especially at large px since valence quarks are distributed 
at relatively large x (partons' fractional momenta of the 
nucleon) while gluons at small x. In other words, high 
Pt protons will have smaller relative contribution from 
gluon jets than anti-protons. If gluon jets lose more en- 
ergy than quark jets as we have assumed in this paper, 
one should then have different suppression factors for 
proton and anti-proton. Such flavor dependence should 
be most evident for heavy particles like nucleons and 
lambdas whose fragmentation functions from a valence 
quark are significantly harder (i.e., falls off more slowly 
at large z) and are very different from gluons and sea 
quarks. For light mesons like pions, the valence quark 
fragmentation functions arc softer and are not much dif- 
ferent from gluons and sea quarks. One then will not 
see much difference between the suppression factors for 
TT"*" and 7r~ even though gluons and quarks have different 
energy loss. 
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FIG. 6. The ratio of 7r~ to 7r+ spectra as functions of 
Pt in pp (dashed) , central Au + Au colhsions at ^/s = 200 
AGeV without energy loss (dot-dashed) and with energy loss 
of dEq/dx = 1 GeV/fm (The mean- free-path Ag = 1 fm). 



Before we discuss the suppression factors, let us look 
at the flavor dependence of the spectra first. Plotted in 
Fig. ^ are 7r~/7r+ ratios as functions of px in pp, cen- 
tral Au + Au collisions with and without energy loss at 
a/s = 200 AGeV. Because gluon-quark scattering domi- 
nates in this Pt region at the RHIC energy and there 
are twice as much valence w-quarks than d-quarks in 
pp collisions, this ratio decrease with pT (dashed line) 
and should saturate at about 0.5 (valence d to u-quark 
ratio in a proton) at very high pT where only valence 
quarks contribute to pion production. At low pT where 
contributions from sea quarks and giuons become more 
important the ratio is then close to one. This is a clear 
prediction of QCD parton model and has been verified 
by experiments some years ago p6[ . In Au + Au colli- 
sions, however, there are slightly more valence cZ-quarks 
than w-quarks since the nuclei are slightly neutron rich. 
As we see in the figure, the 7r~/7r+ ratio (dot-dashed 
line) then increases with pT and approaches to a value 
of about 1.14 which is the valence d to u-quark ratio in 
a Au nucleus. The reason why the ratio is different from 
the limit of d/u ratio is because of finite contributions 
from sea quarks and giuons. If giuons lose more energy 
than quarks, the contribution to highpT pion production 
from giuons will be reduced relative to quarks. There- 
fore, the 7r~/7r+ ratio will be higher than the case of 
no difference in energy loss between quarks and giuons 
(or no energy loss) or become closer in value to the d/u 
ratio, as we see in the figure (solid line). However, the 
change due to the parton energy loss is very small be- 
cause the contributions to pion production from giuons 
is relatively much smaller than quarks. 

The situation for protons and anti-protons is differ- 



ent. From parton distributions in a proton we know that 
gluon to quark density ratio fg/pix,Q^)/fq/p{x,Q^) de- 
creases with X = 2Et/\/^, where Et is the transverse 
momentum of the produced jet. Consequently the ratio 
of gluon to quark jet production cross section always de- 
crease with Et- Since most of anti-protons come from 
giuons while protons come from both valence quark and 
gluon fragmentation, the ratio of anti-proton to proton 
production cross section should also decreases with their 
Pt, as our calculation shows in Fig. ^ for pp collisions 
(dashed line) at ^/s = 200 GeV. At small pt gluon and 
quark jet cross sections become comparable, so the ratio 
p/p should increase. But it will always be smaller than 
1 because there will always be more proton than anti- 
proton in nucleon or nuclear collisions due to baryon 
number conservation (and finite net baryon production 
in the central region even from perturbative QCD calcu- 
lation) . 
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FIG. 7. The ratio of p to p (upper panel) and A to A 
(lower panel) spectra as functions of pr in pp (dashed), cen- 
tral Au+Au collisions at -/s — 200 AGeV without energy loss 
(dot-dashed) and with energy loss of dEq/dx — 1 GeV/fm 
(solid) (The mean- free-path A, = 1 fm). Giuons are assumed 
to loss as twice much energy as quarks. The arrows indicate 
the ratio at low pr < 1 GeV/c from HIJING/BJ estimate 
(with baryon junction model of baryon stopping). 

The dot-dashed line in Fig. ^ is the p/p ratio in cen- 
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tral Au + Au collisions without parton energy loss at the 
RHIC energy. Since Au nuclei are slightly neutron rich, 
one should have less proton production per nucleon from 
valence quark fragmentation than pp collisions. Since 
gluon jet production does not change from pp to Au-\-Au, 
the ratio p/p in Au + Au (without energy loss) is then 
a little larger than in pp collisions. If there is parton 
energy loss and gluons lose more energy than quarks, 
then as we have argued that p/p ratio should become 
smaller than without energy loss (or gluons and quarks 
have the same energy loss), as shown in the figure as the 
solid line. The result and argument is the same for A/A 
ratio as also shown in the lower panel of Fig. 0. To fur- 
ther illustrate this point, we plot in Fig. js] the particle 
suppression factors for proton, anti-proton, lambda and 
anti-lambda as functions oipT- Because of the increased 
energy loss for gluons over quarks, the suppression fac- 
tors for anti-protons and anti-lambdas is then smaller 
than protons and lambdas. This could be easily verified 
if one can identify these particles at high pT in experi- 
ments. 
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FIG. 8. Particle suppression factors for p, A (solid) and 
p, A (dashed) as functions of pr in central Au -\- Au collisions 
&t y/s = 200 AGeV with energy loss of dEq/dx = 1 GeV/fm 
and mean-free-path A, = 1 fm. Gluons are assumed to loss 
as twice much energy as quarks. The arrows indicate the 
suppression factors at low pr < 1 GeV/c if the soft particle 
production is assumed to have an A^'^ scaling. 



In the calculation of high pT baryon spectra in Figs. ^ 
and H one has to use parametrized fragmentation func- 
tions for baryons similarly to these of mesons | p2| , p^ . 
Though baryon production from jet fragmentation in 
e+e" and e~p collisions has been studied |^^, we could 
not find any parametrized form including the evolu- 
tion. Since Lund model has been proven to reproduce 
the experimental data well, here we use the baryon frag- 
mentation functions parametrized from the Monte Carlo 
simulation of Lund model (JETSET) The parame- 
terizations are given in the Appendix. 

In the calculation of the particle ratio in Fig. ^ and 
and particle suppression factors in Fig. ^ we included 
only contributions from perturbative hard processes. As 
we have discussed before there will also be particle pro- 
duction from non-perturbative processes. These soft par- 
ticle production which are dominant at low pt are not 
likely to change the ix^ jix^ ratio much. However, it 
might change the pjp and A/A ratio, as recent heavy- 
ion experiments show that there is significantly more 
baryon stopping than cither pQCD calculation or a sim- 
ple Lund string model of nuclear collisions. There are 
many models of non-perturbative baryon stopping in nu- 
clear collisions . To take into account of this non- 
perturbative baryon stopping, a baryon junction model 
^\ has been implemented into the original HIJING 
model to describe the observed baryon stopping at SPS 
energy [||. This version of HIJING (or HIJING/BJ) 
model at the RHIC energy gives a ratio of pjp = 0.67 
and A/ A — 0.75. These values should serve as an es- 
timate of the particle ratio at small pT < \ GeW /c, as 
indicated by the arrows in Fig. ^. This then gives us an 
upper bound of the uncertainty for the ratio at around 
Pt ~ 2 GeV/c, depending on the interplay between per- 
turbative and non-perturbative contributions. Similarly, 
the suppression factors at low pT also depend on the A- 
scaling of low-p^ particle production. If we assume an 
A^-^ scaling like we did for all charged particles, then the 
suppression factors at low pT < I GeV/c should be as 
indicated by the arrows in Fig. ^ Since contributions 
from the stopped baryon scale like A, the low pT limit 
for baryons will then be smaller than anti-baryons. This 
is an upper bound of the uncertainty one should bear in 
mind at intermediate pr- At large pT these uncertainties 
will become very small. 

Similarly as we have discussed in the previous sections, 
the particle ratio and suppression factors will all depend 
on the other parameters of the energy loss and its energy 
and A dependence. But these will not change the qual- 
itative feature of the flavor dependence of the particle 
suppression due to different energy loss suffered by glu- 
ons and quarks. Because of finite net baryon density in 
the central region, the baryon and anti-baryon absorp- 
tion in the hadronic phase will be slightly different which 
might also give rise to different suppression factors for 
baryons and anti-baryons. Detailed study of this effect is 
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out of the scope of this paper. However, at very large pr, 
the physical baryons might only be formed outside the 
dense region of hadronic matter. Before then, the color 
neutral object might have very small cross section with 
other hadrons which have already been formed. Thus 
the effect of baryon annihilation might be very small at 
large px- Study of the preliminary data for high pT par- 
ticle production in Ph -I- Ph collisions at the SPS energy 
support this scenario p^. 



VI. CONCLUSIONS 

A systematic study of the effects of parton energy loss 
in dense matter on the high pT hadron spectra in high- 
energy heavy-ion collisions has been carried out in this 
paper. We found the hadron spectra at high pT is quite 
sensitive to how the large- E't partons interact with the 
dense medium and lose their energy before they fragment 
into hadrons, leading to the suppression of high-pT par- 
ticles. The suppression factor as a function of pT is also 
sensitive to the energy dependence of the parton energy 
loss. Even though the nonlinear length dependence of 
the energy loss as suggested by a recent theoretical study 
[0 leads to stronger suppression, one cannot unambigu- 
ously determine the nonlinearity by varying the system 
size. We also studied the flavor dependence of the parti- 
cle spectra and the suppression factor and found that it 
is a good probe of the energy loss, especially the differ- 
ence between the energy loss of a gluon and a quark. 

Because of our lack of quantitative understanding of 
energy loss of the produced high Et parton jets inside 
the dense matter in heavy- ion collisions, our phenomeno- 
logical study in this paper can only be qualitative. But 
even such a qualitative study is essential to establish 
whether there is parton energy loss at all in heavy-ion 
collisions and thus whether there is such an initial stage 
in the collisions when the produced dense matter is equi- 
librating. The analysis we proposed in this paper, which 
is also somewhat model dependent, can at least provide 
information about the average total energy loss the par- 
ton could have suffered during its interaction with the 
medium. Anything beyond that will require our knowl- 
edge of the dynamical evolution of the system. Even 
toward such a modest goal, there is still one final hurdle 
to overcome, i.e., final state interactions between lead- 
ing hadrons of a jet and the soft particles in the hadronic 
matter. Such an issue is very important to the determi- 
nation whether the high-p^ particle suppression, if any, 
is indeed caused by parton energy loss in the initial stage 
of a dense partonic matter. Since the formation time 
of a large pt particle is longer than the soft ones, a 
physical large-p^ particle might be formed outside the 
dense region of the hadronic phase. Before then, it is 
in a form a color dipole which might have very small 
interaction cross section with other hadrons. Therefore, 



the hadronic phase of the dense matter might have very 
small effect on the high-pr particle spectra. One could 
address this issue in heavy-ion collisions at the SPS en- 
ergies [|l^ , where one would at least expect that a dense 
hadronic matter has been formed. 



APPENDIX 

In this Appendix we list the baryon fragmentation 
functions parametrized from the results of Lund JET- 
SET Monte Carlo program We simulate the frag- 
mentation of a or a two-gluon system with invariant 
mass W = 2(5, and then parameterize the particle distri- 
butions along one direction of the jet axes as functions of 
z = Eh/Q. We choose the form of the parameterization 
as 

D'liz, Q) = K(Q))A^z"(l - zf{l + z)\ (14) 
if the parton a is a gluon or sea quark, and 

D'liz^Q) = K(0))[A^iz"ni - -^f'i^ + -^V' 

+ N2Z'''{1- zf% (15) 

if the parton is a valence quark of the hadron h. The frag- 
mentation functions are normalize as dzD'^{z,Q) = 
{nh{Q))- For a rough approximation which is enough for 
a qualitative study in this paper, we neglect the change 
of the shape of distributions according to the QCD evolu- 
tion and attribute the energy dependence to the average 
multiplicity {nh{Q)), which are parametrized as. 



{nh{Q)) = a + bs + cs^ 
s = In 



ln(QVA|cD) 



where we choose Qo 



In(g2/AQCD)' 
1 GeV. 



(16) 



1. Gluons: 



a = 0.061, 6 = 0.147, c = 0.155 
iV = 3.814, a = -0.187, /3 = 3.660, 
7 = -2.231 (17) 



a = 0.047, 6 = 0.161, c = 0.133 
iV = 3.814 a = -0.187, /? = 3.660, 
7 = -2.231 (18) 



a = 0.0215, 6 = 0.0454, c = 0.0568 
TV = 3.378, a = -0.166, /3 = 4.394, 
7 = 0.105 (19) 



10 



2. d quarks: 



a = 0.0966, h = 0.0419, c = 0.1045 
N2 = 1.671, a2 = 0.699, (32 = 1-311 
Ni = 0.002, ai = -2.303, /3i = 6.461, 
71 = 20.225 (20) 



a = 0.0392, h = 0.0356, c = 0.0906 

N2 = 1.377, Q!2 = -0.252, /32 = 2.142 

Ni = 0.005, ai = -2.246, /3i = 4.464, 
71 = -2.141 (21) 



£)2 



a 
AT 

7 



a 
N 



0.0098, 6 = 0.0269, c = 0.0250 
0.230, a = -1.027, /3 = 1.962, 
3.037 (22) 



0.0104, b 
0.318, a = 
5.186 



0.0867, 
-0.989, 



c = 



0.0743 
= 4.956, 



(23) 



a = 
N = 

7 = 



0.0124, 6 = 0.0676, c = 0.0760 
0.318, a = -0.989, /? = 4.956, 
5.186 (24) 



a 
N 



0.0033, 6 
0.318, a = 
5.186 



0.0232, c = 0.0265 
-0.989, /3 = 4.956, 



(25) 



3. u quarks: 



By isospin symmetry: Du'"^^ = Z)^^"\ Z)"^"^ = 



4. s quarks: 



D 



A(A) 



a = 0.0706, b = 0.0546, c = 0.0113 
ATs = 11.880, a2 = 2.790, (32 = 1.680 
ATi = 9.55 X 10-^ ai = -3.09, /3i = 8.344, 
71 = 31.74 (26) 



a = 0.0362, b = 0.0228, c = 0.1087 
AT = 0.254, a = -1.0123, /? = 3.506, 
7 = 4.385 (27) 



a = 0.0326, b = 0.0149, c = 0.1060 
N = 0.421, a = -0.867, /? = 3.985, 
7 = 3.577 (28) 



a = 0.0123, 6 = 0.0631, c = 0.0869 
Af = 0.410, a = -0.931, /3 = 5.549, 
7 = 4.807 (29) 



• DP 



a = 0.0135, 6 = 0.0456, c = 0.0935 

iV = 0.410, a = -0.931, /? = 5.549, 
7 = 4.807 (30) 



a = 0.00197, b = 0.0331, c = 0.0174 
TV = 0.238, a = -1.060, /? = 7.141, 
7 = 9.106 (31) 

5. Aiiti-quarks: By symmetry of charge conjugate: 

q q q q 
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